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ABSTRACT 

Type la supernovae are believed to be wiiite dwarfs disrupted by a tiiermonuclear explosion. Here we investigate the scenario in 
which a rather low-mass, carbon-oxygen (C -I- O) white dwarf accumulates helium on its surface in a sufficient amount for igniting a 
detonation in the helium shell before the Chandrasekhar mass is reached. In principle, this can happen on white dwarfs accreting from 
a non-degenerate companion or by merging a C + O white dwarf with a low-mass helium one. In this scenario, the helium detonation 
is thought to trigger a secondary detonation in the C + O core. It is therefore called the "double-detonation sub-Chandrasekhar" 
supernova model. 

By means of a set of numerical simulations, we investigate the robustness of this explosion mechanism for generic 1-Mq models and 
analyze its observable predictions. Also a resolution dependence in numerical simulations is analyzed. 

Hydrodynamic simulations of the double-detonation sub-Chandrasekhar scenario are conducted in two and three spatial dimensions. 
The propagation of thermonuclear detonation fronts, both in helium and in the carbon-oxygen mixture, is computed by means of both 
a level-set function and a simplified description for nuclear reactions. The decision whether a secondary detonation is triggered in the 
white dwarf's core or not is made based on criteria given in the literature. 

In a parameter study involving different initial flame geometries for He-shell masses of 0.2 and O.I Mq (and thus 0.8 and 0.9 Mq of 
C + O), we find that a secondary detonation ignition is a very robust process. Converging shock waves originating from the detonation 
in the He shell generate the conditions for a detonation near the center of the white dwarf in most of the cases considered. Finally, we 
follow the complete evolution of three selected models with 0.2 Mq of He through the C/O-detonation phase and obtain ^*Ni-masses 
of about 0.40 to 0.45 Mq. 

Although we have not done a complete scan of the possible parameter space, our results show that sub-Chandrasekhar models are not 
good candidates for normal or sub-luminous type la supernovae. The chemical composition of the ejecta features significant amounts 
of ""^Ni in the outer layers at high expansion velocities, which is inconsistent with near-maximum spectra. 

Key words, supernovae: general - nuclear reactions, nucleosynthesis, abundances - hydrodynamics - methods: numerical 



1. Introduction 

One of the major uncertainties in modeling type la supernovae 
(SN la) originates from the unknown nature of the progenitor 
systems because neither observations nor theoretical models are 



yet conclusive (Branch et al. 1995[ |Ruiz-Lapuente et al.| | 2000 



|Livio"2000', 'Nomoto et al."2003', 'Han & Podsiadlowski"2004 
Napiwotzki et al. 2005; Stritzinger et al. 2006; Parthasarathy] 
et al.||2007| l. Over the past years, most studies of thermonu 



clear supernova explosions focused on models in which a ther 
monuclear flame is formed by a runaway near the center of a 
white dwarf (WD), composed of carbon and oxygen, once it 
has reached the Chandrasekhar mass (Mch ~ 1.4 Mq) by ac- 
cretion from a non-degenerate companion. Starting out in the 
sub-sonic deflagration mode of flame propagation, these mod- 
els were shown to give rise to energetic explosions either in 
pure turbulence-booste d deflagrations (|Re inecke e t al.||2002b 



Gamezo et al. 2003 ; Ropke & Hillebrandt 2005 , Ropke et all 



2006 ; Schmidt & Niemeyer 2006; Schmi dt et al.f2 006; Ropke 
et al.,|2007a^ or with a delayed triggering of a supersonic det- 



onation phase (Gamezo et al. 2004[ Plewa et aL]|2004[ [Ropb 



|& Niemeyer|2007| ). There seems to be reasonable agreement of 
such models with the gross features of observed SNe la ( Ropke ^ 
|et al.|[20"07a[ [Mazzali et al.||20(j7l ). However, it remains unclear 
whether the Chandrasekhar mass model can account for the full 



range of SN la observations. In particular, the mechanism of the 
sub-luminous events (SN 1991bg being a prototypical example) 
remains a puzzle ( |Mazzah et al. 2007). St ritzinger et al .1 ( 120061 1 
even claim typical ejecta masses below Mqu for a sample of well 
observed SNe la, with the trend that the low-luminosity explo- 
sions eject less mass. 

Therefore, a long standing question is whether other pro- 
genitor channels contribute to the observed SN la sample, and 
two alternatives have been suggested: the WD merger or dou- 
ble degenerate scenario (in contrast to single degenerate mod- 
els, which consist of binaries with only one WD) and the sub- 
Chandrasekhar model. In the present work we explore the latter. 

The basic idea of the sub-Chandrasekhar model is that if the 
accretion rate onto a WD is lower than about 1 - 4 ■ 10"'^ Mq yr"' 
the accreted He (or the H processed to He) is not fused steadily 
into C and O. Instead, after reaching a critical amount of He 
at relatively low densities, the He shell becomes unstable and 
detonates (cf. [Woosley et al.|1986| and references therein). The 
detonation ignites most Ukely close to the bottom of the He 
shell, produces almost pure ^^Ni, and can occur long before the 
WD reaches the Chandrasekhar limit. A second detonation may 
be triggered spontaneously when the He shell detonation shock 
wave hits the C -i- O core or with some delay after the shock has 
converged near the center of the WD (double detonation mod- 
els). This way, in the sub-Chandrasekhar mass models the condi- 



2 



M. Fink et al.: Double-detonation superaovae of sub-Chandrasekhar mass white dwarfs 



tions for thermonuclear runaway are caused by the compression 
in the shock wave and not by the high degree of degeneracy as 
in the Chandrasekhar mass models. 

Delayed double detonations have been studied extensively 
before. [Woosley & Weaver] ( [T994| , and 
one-dimensional (ID) simulations and 



Livne ( 1990 1 carried out 
Livne & Glasnerl ( fT990l 



11991 ) considered two-dimensional (2D) setups. In ID the He 
detonation ignites synchronously in a layer close to the bottom 
of the He shell and due to the spherical symmetry constraint the 
shock converges perfectly in the center. However, this is not a 
very realistic model. According to Livne ( 1990) an ignition most 
likely happens in a single point leading to an off-center conver- 
gence of oblique shock waves in the core. 

Several other simulations predict successful directly ignited 
double detonations in two dimensions ( Livne| 1 997) |Arnett| 1 997 [ 



[Wiggins & Falle | 1997 |[ Wiggins et al.|199 8|l. Th e sm oothed parti 
cle hydrodynamics simulations by Benz (1997 1 and Garcia-Senz 
let al.| ( fT999) were carried out in 3D. Livne & Glasner] ( |1990l l 
first reported an increased probability of the direct core ignition 
if the He detonation does not happen directly at the core-shell 
interface but at a certain distance above it. This way, the pressure 
jump can grow large enough before hitting the core. 

In this work, the results of a (restricted) parameter study are 
presented investigating the possibility of triggering the second 
detonation in two and three dimensions. In all our models the 
WD has a total mass of 1 Mq, but the (C + 0)-core and the He- 
shell masses differ. We compute sequences of models with very 
different initial flame geometries in order to test the robustness of 
this explosion mechanism. Since this latter question is the main 
focus of our paper, most of the simulations are stopped once 
the conditions for a detonation are matched. However, for a few 
successful cases the energetics and nucleosynthesis of complete 
double detonations are computed. In the following chapter the 
model details will be described. Section [3] shows the simulation 
results, which are summarized and discussed in the last part of 
the paper. 



2. Explosion model 

2.1. General setup 

The numerical scheme used in our simulations is a modified ver- 
sion of the code described b y |Reinecke et al.|(|1999a| 



Ropke & HiUebrandt] poos) [Golombek & Niemeyer 



2002a I 



(20051, 



and Ropke & Niemeyer ( 2007| l. The main difference is the burn 
ing physics which accounts for He detonations and the propaga- 
tion of two detonation waves. 

Our hydrodynamic simulations were carried out on a uni- 
form Eulerian grid with a finite volume scheme based on the 
PROMETHEUS implementation ( [FryxeU et aL][T989| of the 
"piecewise parabolic method" (PPM) (Colella & Woodward 
1 1984) 1. Thus, discontinuities like shocks were automatically and 
accurately captured but, nonetheless, smeared out over several 
grid cells. A potential concern is the overheating effect ob- 
served in Godunov-type schemes when shocks hit walls or col- 

{T996|, 



lide. Following the suggestion of Donat & Marquina 



an approximate Riemann solver (the so-called Marquina solver) 
was implemented but this did not lead to noticeable differences. 
Thus, all simulations presented in the following are based on the 
exact Riemann solver of |Colella & Glaz| ( |1985) l. 

While the parameter study testing the possibility of igniting 
a detonation in the C -i- O core was carried out on a static grid, 
the complete double detonation models presented at the end of 
Sect.[3)employed a co-expanding grid as in |R6pke & Hillebrandt| 



( |2005) l in order to account for the expansion of the WD. The 
coordinates were (r, z)-cylindrical in 2D and cartesian in 3D. 

The equation of state we used includes contributions from an 
arbitrarily degenerate and arbitrarily relativistic electron gas, a 
photon gas, electron-positron pairs and nuclei with a Maxwell- 
Boltzmann distribution. The initial WD was assumed to be in 
hydrostatic equilibrium. Therefore its structure is uniquely de- 
termined once the density at the center and at the core-shell in- 
terface and the temperature profile are fixed. For simplicity a 
constant temperature of 5 ■ 10^ K was assumed throughout the 
star Of course, this is an unrealistic choice, but because of the 
high degree of degeneracy the exact temperature has almost no 
effect on the state of the matter. Another simplifying assumption 
is the approximation of the gravitational potential in spherical 
symmetry. Only the monopole moment of the density distribu- 
tion averaged over all angles was calculated at every time step. 

2.2. Nuclear reactions 

As we are mostly interested in the hydrodynamic evolution 
of the explosion, it is appropriate to use a strongly simplified 
scheme for nuclear reactions as in |Reinecke et al.| ( |2002a[ ). On 
our numerical grid the structure of the flame cannot be resolved. 
Therefore in our simulations the detonation flame is approxi- 
mated as a discontinuity propagating at an appropriate speed (see 
Sect. 2.3 1 with the fast reactions happening spontaneously at its 



passage and only the density dependent products, but no reac- 
tion rates being calculated. In order to approximate the energy 
release in the fast reactions and the correct molecular weight of 
fuel and ashes, it is sufficient to consider only one representative 
species for each group of elements involved which has a bind- 
ing energy and mass typical for the group. Our scheme includes 
four species, namely "^He, a mixture of '^C and "'O with equal 
mass fractions, ^'^Si as a representative for the intermediate mass 
elements, and ^^Ni as an iron group element. The following re- 
actions are included in the model, depending on the density at 
the front: 

- The ^He detonation in the shell produces pure ^''Ni. Flame 
extinction is assumed to occur atp = 10^ gem"-'. 

- A detonation in the '^C/"'0 core produces ^^Ni in nuclear 
statistic equilibrium (NSE) (represented by a temperature 
and density dependent mixture of ^''Ni and '^He) above p = 
3-10^ gcm"^* and -'^Si below^At densities smaller than 
lO*" g cm"-' burning is assumed to cease. 

Since the transition density to NSE affects the ''^Ni-mass and 
thus the total explosion energy, its exact value should be investi- 
gated further in the future. 

In order to distinguish the "^He in the shell from that produced 
in NSE with ^''Ni both are treated as different species in the code. 
With decreasing densities the "^He in NSE finally recombines to 
^^Ni. This is the only slow reaction included in our simulations. 

The propagation of the detonation flame is modeled with the 
"level set technique " (fOsher & S ethian 1988) in the implemen- 
tation of |Golombek & Niemeyer| ( |2005), which in turn is based 



on the deflagration flame model of Reinecke et al. ( 1999b) ). The 



zero level set of a signed distance function G represents the flame 



According to 



Imshennik & Klioklilov 



1984 1, at densities lower 



than about 10' gem"" the timescale for NSE formation becomes so 
large that the flame width would be bigger than the extent of the whole 
core. Since it can be assumed that NSE is no longer reached even before 
the flame width gets close to the extension of the core, 3 • 10' gcm"^ 
instead of 10' g cm"'' is used as the transition density. 
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front. Information on whether matter has already been burnt is 
encoded in its sign. It is negative in the fuel and positive in the 
ashes. Most important, the location of the flame can be calcu- 
lated at every time step simply by finding the roots of G. The 
propagation of the zero level set of G consists of two parts: the 
advection with the flow and self-propagation by burning, respec- 
tively. By means of an operator splitting, the advection of G can 
be treated in the same way as that of a passive scalar in the hydro 
dynamic scheme (Mul der et al.|1 992, Rei necke et al.|1999a|b 
Burning is accounted for by ( Reinecke|2001 



IdG 



=s\VG\. 

\ /burn 

Thus the velocity s of the flame relative to the flux has to be de- 
termined at every time step from the quantities on the grid. In the 
state-of-the-art deflagration codes this is done by an appropriate 
sub-grid turbulence model. For detonations, as occurring in our 
models, turbulence effects can be neglected and the propagation 
speed s can be determined easily as a function of the density p 
only, as will be described in the next section. 

The changes in the mass fractions are calculated depending 
on p, when the flame front passes a cell. Here and also in the 
case of the velocity s we have the problem that in general the 
state in such a cell is a mixture of a burnt and an unburnt matter 
So we have e.g|^ 

p^api+{\- a)p2, 

where a is the volume fraction of the cell occupied by the un- 
burnt state. With the so called "complete coupling" scheme the 
two partial states could be reconstructed from the mixed state 
( Smiljanovski et al. 1997 1. This enables flux splitting that avoids 
numerical mixing of fuel and ashes. While such a scheme has 
been implemented for 2D small-scale deflagration flame simula- 
tions in WD matter ( Ropke et al.|2003 2004a|b[ l, numerical ob- 
stacles prevent its application in full star simulations up to now. 
Here, we apply the "passive implementation" of Reinecke et al.| 
( | 199 9b ) and use the mixed state to advect G and also to calculate 
the burning velocity s and to determine the reaction products]^ 
The scheme described here is not in general appropriate to 
simulate detonations as, e.g., the burning speed s depends on the 
detailed geometry of the detonation front. The burning speed s is 
in general a function of |, where 6 is the width of the detonation 
front and R is its curvature radius. In our case, however, the front 
can be approximated as thin and planar in the density range that 
is relevant for the explosion dynamics. 



2.3. Detonation physics 

A detonation wave consists of a shock wave that compresses 
and heats the matter behind it so strongly that it is ignited. The 
shock wave, in turn, is maintained by the energy released due 
to burning. Depending on the properties of the reactions taking 
place, there exist different types of detonations. A very impor- 
tant case is the Chapman-Jouguet detonation, which propagates 



into the fuel at the lowest possible speed. According to Landau 
|& Lifschitz| ( |199l| a detonation that is created spontaneously by 
the combustion itself is a Chapman-Jouguet detonation in many 



^ Here and in the following indices 1 and 2 refer to quantities on the 
unburnt and the burnt side of the front, respectively. 

3 Correctly, one would have to use burning and flow velocities on 
either the burnt or the unburnt side of the flame to update G due to 
burning and advection. 



Table 1. Critical masses for detonations at given density and 
temperature according to Niemeyer & Woosley ( 1997| l. 



P [gem 3] 


[10'-* K] 


M[g] 


R[m\ 


2-10^ 


2.8 


2.0- 10'" 


0.7 


10« 


3.2 


2.0-10'" 


2 


3 • iC 


3.2 


2.0- 10" 


50 



Table 2. Critical temperatures for detonations at given density 

( 2007b[ l. Detonation is im- 



and mass according to 'Ropke et al. 



possible forp < 1 ■ 10 gem 



P [gem 3] 


[10'' K] 


M[g] 


R [km] 


10' 


2.8 


2.5 • 10-3 


2 


10' 


2.0 


2.0 • 10-" 


8 


10' 


1.9 


1.5- 10" 


30 


3- 10" 


2.3 


2.0 • lO-** 


120 



cases. Therefore, in our simulations we assumed this type of det- 
onation for the He burning. In the Chapman-Jouguet case there 
exists a simple law for the velocity of the burning front relative 
to the ashes, which is just equal to the sound speed there: 



«2,He = 



\ op2 It 



72 



Here y = ^ is the ratio of the specific heat at constant pressure 
and the specific heat at constant volume. 

According to 'Sharpe ( 1999| l, self supporting detonations in 
degenerate carbon/oxygen mixtures at densities between 2-10' 
and 1 ■ 10^ g cm""* are of pathological type. These arise in general 
from a non-monotonic heat release in the reactions. For them the 
front velocity cannot be determined independent of the reaction 
rates. Therefore, for the C/O detonations in this work, we used 
the values calculated by Sharpe ( ^1999^ in numerical simulations 
as a function of density: 



■Sl.C/O - /Sharpe(p)|^ 



2.4. Ignition conditions in tine C + O core 

As our code lacks a detailed nuclear network and reaction rates 
are not calculated, the onset of a C/O detonation can not be fol- 
lowed. Apart from that, in our full star simulations it would not 
be feasible to resolve the scales of the detonator (where the shock 
that supports the detonation is built up). Therefore, the decision 
on whether a detonation is triggered is made upon critical con- 
ditions given by "Nieme yer & Woosley] ( |1997[ ) and [Ropke et al.| 
(2007b ) (see Tables [T] and |2|). They are based on the follow- 
ing arguments. If there is a sufficiently large region of super- 
sonic spontaneous burning, a strong pressure wave must form 
and will synchronize with the burning wave if the compression 
and heating is strong enough to burn the material in the sound 
crossing time of the region under consideration ( |Blinnikov &| 
|Khokhlov|[T986l |Woosley||T990t |Khokhlov|[T99T] l. Therefore a 
detonation wave is bom if that region is big enough or, equiva- 
lently, provided that it consists of enough mass at a given density. 



■* We neglect here that over-driven detonations powered by shock 
waves entering the core due to the He shell detonation could also play a 
role. 
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|Niemeyer & Woosley| ( |1997l l empirically determined those crit- 
ical masses from hydrodynamic simulations for different densi- 
ties and for a given lineaflradial temperature gradient across that 
mass Pj Ropke et al.| (|2007b ) used the same code and setup, but 
determined the critical temperature at a given mass for different 
densities. They focussed on the low end of possible densities. 



3. Simulations 

With the numerical scheme described above, several simulations 



were performed. In Sect. 3.1 the results of a parameter study 
on the influence of the initial flame geometry for two different 
core and He-shell masses are presented, followed by a resolution 
study in Sect. 



3.2 



and, finally, the determination of explosion 
energies and nucleosynthesis yields in three successful double 
detonation simulations. We stress again that the main intension 
of the parameter study is not to scan the entire parameter space 
of possible WD properties but rather to investigate the conditions 
which may lead to a double detonation. 



3. 1. Initial flame geometry 

The accretion of matter from the companion star onto the WD is 
a very complex process, and the formation and ignition of a He 
shell is still afflicted with large uncertainties. As many different 
initial flame geometries in the shell seem possible, their influ- 
ence on the possibility and the properties of a double detonation 
is explored first. 

The different initial flame models, which were treated in this 
study, are shown in Fig. [T] and the nomenclature of the simu- 
lations is as follows: The first letter in each model's name de- 
notes the geometry, and the following number specifies it fur- 
ther. "s4.10" means an ignition in a whole layer of the He shell 
at a radius of 4100 km. "b«" stands for an ignition in n separate 
"bubbles" (which are, of course, tori in 2D if located off the z- 
axis). Two special cases of the latter type, for n = 1, are "z5.00" 
and "rS.OO", where the initial "bubble" is located at the z- or the 
r-axis, respectively, at a distance of 5000 km from the origin. 
Finally the capital letters "A" or "B", which follow the above 
names, indicate the initial masses which were used: 



- Model A: 

- Model B: 



Mc/o = 0.8 Mo 

= 0.9 Mr, 



M. 



c/o 



MHe = 0.2 Mo, 
MHe=0.1Mo. 



The results of the study are given in Table |3] There the di- 
mensionality and the number of grid cells along the side of one 
quadrant are added to the model names. Throughout the study 
a resolution of 256 cells was used. "2dq" means that only the 
first quadrant has been simulated assuming equatorial symme- 
try, whereas "2d" denotes full two dimensional simulations of 
the whole z-range, capturing the entire WD star 

The most symmetric case with the simplest hydrodynamic 
flow pattern is the ignition in a whole layer of the He shell. In 
model s4. 10A_2dq_256 starting with the initial flame geometry 
shown at the top left of Fig.[T] the inner and the outer spherical 
detonation front start to propagate inwards and outwards, respec- 
tively. While the combustion stops close to the border of the WD 
due to the low densities, it ceases at the core-shell interface be- 
cause of the depletion of the He fuel. To investigate the chances 



^ with respect to inner mass 

^ In Table 111 only the results for a mixture of equal mass fractions of 
C and O are shown. In their work |Niemeyer & Woosle3^ ( |1997^ addi- 
tionally investigated the dependence on the mass fraction. 



Table 3. Results of the parameter study about the initial flame 
geometry. The last column states if the ignition conditions of a 
core detonation have been reached. 



Model 


f/s 


-^9. max 


Pi, max 


Core ignition 


s3.60A_2dq_256 


0.627 


9.09 


6.14 


yes 


s4.10A_2dqJ56 


0.575 


10.2 


131 


yes 


s5.00A_2dq_256 


0.609 


12.1 


9.05 


yes 


r4.24A_2dq_256 


0.699 


5.62 


1.54 


yes 


r5.00A_2dq_256 


0.718 


6.00 


1.69 


yes 


z4.24A_2dqJ56 


0.989 


3.27 


1.03 


no 


z5.00A_2dq_256 


1.03 


3.20 


0.981 


no 


blA_2dq_256 


0.767 


4.43 


1.32 


yes 


b2A_2dq_256 


0.678 


8.61 


4.85 


yes 


b3AJdq_256 


0.674 


7.94 


4.57 


yes 


b5iA_2dq_256 


0.629 


9.53 


6.68 


yes 


b5iiA_2dq_256 


0.648 


9.00 


6.13 


yes 


b5iiiA_2dq_256 


0.648 


9.03 


6.37 


yes 


z4.24A_2d_256 


1.08 


8.53 


4.52 


yes 


s4.30B_2dq_256 


0.791 


6.82 


4.69 


yes 


s4.80B_2dq_256 


0.745 


8.17 


5.99 


yes 


r4.50B_2dq_256 


0.892 


3.96 


1.15 


no 


r5.00B_2dq_256 


0.885 


4.49 


1.23 


yes 


z4.50B_2dq_256 


1.21 


2.42 


0.808 


no 


z5.00B_2dq_256 


1.24 


2.39 


0.768 


no 


b2B_2dq_256 


0.855 


6.72 


4.34 


yes 


b5iB_2dq_256 


0.806 


7.63 


5.47 


yes 


z4.50B_2d_256 


1.31 


6.50 


3.38 


yes 



of a core ignition, only the hydrodynamic evolution was consid- 
ered afterwards and all the reactions taking place were stopped. 
The shock wave that was formerly driving the inner detonation 
front continues to propagate inwards until it converges at the 
center. There (at t - 0.575 s) a strong shock collision takes place 
that causes a maximum temperature of 10.2 ■ 10'^ K and densities 
up to 7.37 ■ 10^ gcm"^ 

According to Table [T] these maximum values suffice for a 
C/O ignition, if there is at least a mass of about 2.0 ■ 10'^ g 
with the same temperature and density or, more accurately, with 
a sufficiently small temperature decrease from the above men- 
tioned value. A mass of 2.0 ■ 10'^ g at this density corresponds 
to a spherical volume with a radius of the order of 1 km. This 
length was not resolved in this simulation, where one cell has 
a diameter of 23.2 km. But as the high temperatures and den- 
sities in our simulations are not confined to one single cell, but 
distributed over several (which means that the region of interest 
is resolved), the ignition conditions seem to be clearly fulfilled. 
Consequently, a C/O detonation will occur at the center of the 
core. There is, however, the possibility that a detonation of the 
C/O material is triggered already when the shock wave enters at 
the edge of the core. About 0.1 s after the start of the simula- 
tion, temperatures of about 1 to 2 ■ 10^ K and densities of about 
10^ gcm"^ are reached. According to Table [2j this is within the 
range of conditions where a detonation might occur But as the 
temperature gradient might possibly be too large across a radius 
of about 30 km, this question remains open here. As was already 
mentioned, other aufliors (|Livne|I997| 'Arnett'1 997[[Benz|I997 



Wiggins & Fane| |I997| [Wiggins et al. 1998; G arcia-Senz et al. 
1999j) found successful direct ignitions and an increasing prob- 
ability for the latter with increasing distance of the He ignition 
point from the C/O core (cf. Livne & Glasner|I990| l. 

Our simulations confirm the earlier results. Helium shell det- 
onations were ignited at smaller (s3.60A_2dq_256) and larger 
(s5.00A^dqJ256) radii than in the case described above. As was 
expected, the shock wave propagating through the core is weaker 
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He 


S4.10A 


He 


Z4.24A 


He 


r4.24A 


C + O 


H , 1 f-^ , ^^ 
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— 1 — t — t — h- 
He 


blA 


— 1 — 1 — 1 — 1 — 
He 
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— 1 — t — t — h- 
He 
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• 




■ 
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• 


C + O 
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C + O 
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C + O 
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H , 1 1 , , ^^ 


— 1 — < — < — h- 
He 


b5iA 
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• 
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r[1000kmj 



2 4 
/■ [1000 km) 
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Fig. 1. Initial flame models. The He detonation flame fronts start at the borders of the black regions that already contain hot ashes at 
the beginning. 



in the case of smaller ignition radii and stronger in the case of 
shocks starting further out. Accordingly the maximum tempera- 
tures and densities that are reached in the center are lower in the 
former and higher in the latter case (see Table|3]l. 

Due to the high degree of symmetry, the models with shell 
ignition are the ones that can trigger a second detonation most 
easily. To test multi-dimensional configurations, the complete ig- 
nition layer was replaced by single spots (which effectively are 
tori in 2D) distributed over the He shell. Although this is proba- 
bly not a realistic scenario, it is considered here as a numerical 
experiment. As far as their maximum densities and temperatures 
are concerned, the models with five initially ignited bubbles in 
the first quadrant differ only slightly from the complete shell ig- 
nition model (see Table|3]l. The values of Tmax and p^ax are only 
slightly lower than for the ignition of a shell, starting at the same 
radius, and they still match the conditions for a detonation. If 
the number of initially burning bubbles is reduced to three, peak 
temperature and density decrease further 

The evolution of simulation b3A_2dq_256 is shown in Fig.|2] 
Due to the coalescence of the parts of the detonation fronts mov- 
ing inwards, an almost spherically symmetric shock front is en- 
tering the core. This resemblance of the dominant hydrodynamic 
processes is the reason for the similarity of the results given in 
Tabled 

In between the bubbles interactions of the detonation shocks 
take place. The collision process of two radial shocks can be 
described as a reflection at the plane of symmetry between the 
waves (see e.g. [Courant & Friedrichs||1948) l. Generally, in this 



frame the incidence angle may differ from the angle of reflection. 
This can be seen best at f = 0.15 s in Fig.|2] The reflected waves 
are again reflected at the waves coming from the opposite side of 
the initial bubble. The repetition of this process should also have 
some influence on the geometry of the innermost shock front and 
thus on the conditions given in Table [3] 

This can be seen in model b2A_2dq_256, which is simi- 
lar to model b3AJ2dq^56, but with the bubble on the z-axis 
removed. It starts with only two initially ignited bubbles but 
reaches greater maximum densities and temperatures in the core. 
The reason may be a focusing effect on the z-axis in 2D cylindri- 
cal symmetry and the bigger density and pressure jump that can 
build up until the burning reaches the axis. 

An analogous focusing effect in the plane of equatorial sym- 
metry (the xy-plane) is not observed. On the contrary, remov- 
ing the initial burning bubble from the r-axis (i.e. going from 
b2A_2dq^56 to blA^dq_256) leads to a strong decrease of the 
maximum values of density and temperature. 

Considering models with only a single initially burning bub- 
ble, the exact geometry of the appearing shock waves is of great 
importance. After the He detonation has reached the border of 
the core, it spreads around it in a wave propagating almost per- 
pendicular to the core-shell interface. The detailed hydrody- 
namic structure has been analyzed in |Livne & Glasner ( 1990[ l, 
from which Fig. [3] is taken. An oblique shock wave that propa- 
gates into the core (tilted by an angle 6 with respect to the inter- 
face) is associated with the detonation wave moving along the 
surface of the core. After the shock passage the core-shell inter- 
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Fig. 2. Time evolution of the model b3A^dq^56. Here and in the following figures the density is color coded and lines of constant 
pressure are drawn in blue. In cyan the solid line is the location of the He detonation flame (the zero level of G) and the dashed hnes 
are at the border of the He shell. 

applying the Rankine-Hugoniot conditions and depend on the 
strength of the detonation shock. 

This structure can clearly be seen at f = 0.40 s in model 
r4.24A_2dq_256, which is the detonation of a torus in the z = 
plane (see Fig. |4]). Starting from the "equatorial" plane the det- 
onation wave propagates around the core until it reaches the 
"polar" z-axis. There a collision of the waves coming from all 
azimuthal directions takes place. Equivalently, each wave is re- 
flected at the wave coming from the opposite side of the star. 
This causes the shock wave visible in the last two snapshots of 
Fig.|4]to propagate to the lower right. In the following the oval 
inner shock front moves towards the center. At f = 0.70 s, even 
before the center is reached, the maximum density and tempera- 
ture values are attained on the z-axis at about 850 km away from 
the center (fourth plot in Fig. |4]). The density and temperature 
values of 1 .54- 10^ g cm"-* and 5.62 ■ 10^ K are significantly lower 
than those in the spherically symmetric models, but according to 
Table [T] they still suffice for the ignition of an off-center core 
detonation. 

In the following simulations, the symmetry of the first deto- 
nation was further reduced by departing from equatorial symme- 
try and simulating a whole rotationally symmetric star, i.e. now 
the negative z-range was considered too. Fig. |5] shows the exam- 
ple simulation z4.24A_2d^56. Here the He detonation starts in a 
point on the z-axis and then spreads around the whole WD until 
it reaches the negative z-axis (see Fig.|5]). In this sense, the setup 
is maximally asymmetric. The processes now taking place are 
analogous to the simulation r4.24AJ2dq^56 described above. 
But here the maximum density and temperature are reached 



HE , (rarefaction) 




Fig. 3. Structure of a detonation propagating at the core-shell 
interface in the rest frame of the detonation front from lLivne &I 
[Glasner ( 1990). (Here the core consists of carbon only, not of 
C/O like in our case.) 



face is also tilted inwards by an angle y6. Thus in-flowing matter 
follows the oblique shock. The detonation wave is followed by 
a rarefaction wave. The angles B and 6 can be determined by 
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t = 0.70 s 



r [1000 km] 

Fig. 4. Time evolution of the model r4.24A^dq_256. 



)■ [1000 km] 



when the oval inner shock front converges almost perfectly in 
a "point" on the negative z-axis a.tz~ -1500 km. Thus at a time 
of 1.08 s, = 4.52 ■ 10** gcm-3 and r^^x = 8.53 ■ lO** K 
are reached. These maximum values are comparable with the 
spherically symmetric case discussed at the beginning of this 
section. The reason for this is the high degree of symmetry of 
the colliding shock fronts. Astoundingly those maximum values 
are reached comparatively far away from the center in a region 
of relatively low initial density. 

The initial condition considered in the last paragraph seems 
more realistic than the more symmetric cases. Thus, this model 
was studied also in simulations with higher numerical resolution. 
Here, a resolution dependence of the peak temperatures and den- 
sities became apparent - an effect that is discussed in Sect. 3.2 
below. 

Finally, a series of models with a different initial mass dis- 
tribution was investigated. Sequence "B" corresponds to mod- 
els with 0.9 Mq of C/O and 0. 1 Mq of He and the parameters 
given in Table [3] They support the generality of the results dis- 
cussed above. The dependence on the initial flame geometry is 
analogous to case "A", but maximum densities and temperatures 
achieved are generally lower This is consistent with the fact that 
the shock waves originating from the He shell have to cross a 
larger pressure gradient before they reach the center in model 
"B"|^ Still, most of the models reach the conditions for core ig- 
nition described in Sect. |2.4| (the "zr" models did not ignite with 
mass distribution 'A" either). 



' Geometrical amplification efi'ects do not play a role at such big 
radii. 



3.2. A resolution study 

The numerical resolution of the parameter study presented above 
was relatively low with a spatial resolution of 23.2 km only. 
Hence, for one of the most interesting models, z4.24A^d_256, 
we performed a resolution study. Globally, we found a good 
agreement of the temporal evolution of the state variables. But 
the maximum densities and temperatures at the central shocks 
collision showed a significant dependence on the numerical res- 
olution of the simulation (see Fig. |6]). With decreasing cell size 
the maximum values increase almost exponentially. This trend 
continues down to the smallest tested cell width (7.4 km) and 
no indication for numerical convergence is found. Given that 
the shock collision is an almost perfect spherically symmetric 
problem, we reduced to one dimension in order to facilitate even 
higher resolutions. The state variables on the r-axis of simula- 
tion s4.10A_2dq_256 at f = 0.25 s linearly interpolated to higher 
resolutions were taken as initial condition. The results of the cor- 
responding one-dimensional runs with increasing resolution are 
presented in Fig. |7] It shows an analogous exponential increase 
of the temperature and density maxima down to a cell width of 
0.72 km. Below that value, however, the temperatures exceeded 
the limits of our equation of state. Here, more micro-physical 
processes would have to be taken into account and could restrict 
further growth of density and temperature. 

But even without this restriction, from analytical consid- 
erations it is clear that the density growth ratio must be lim- 
ited, although it can get very big. In the idealized case of in- 
finitely strong shock reflection at the origin in spherical geom- 
etry (and constant initial conditions of po(r) - 1, uo{r) - -1, 



M. Fink et al.: Double-detonation superaovae of sub-Chandrasekhar mass white dwarfs 




1 2 3 4 5 
r [1000 km] 



3 4 5 
)• [1000 km] 



Fig. 5. Time evolution of the model z4.24A_2d_256. 



and Pair) - 0), its value is after the reflected shock 

(cf. |Glaister|l988j lp] Another problem that can be solved semi- 
analytically is that of a spherical shock wave that starts from in- 
finity, converges to the origin, and is reflected there (cf. G uderley| 
[19421 [Landau & Lifschitz J"WT] [Ponchaut et al. 2006 ). While 
the shock is moving inwards it grows in strength because of 
its surface decrease. Close enough to the center for the strong 
shock approximation to hold, it can be shown that p ~ , 
p = const, and \u\ - with n being a constant between 
and 1 that has to be numerically determined and depends only 
on the ratio of specific heats y. According tojLandau & Lifschitzj 



1991| l, at the time of perfect focusing the density has increased 
by a factor of 20. 1 for y - \ and the reflection is accompanied 



^ For an ultra-relativistic degenerate gas like in a WD's center y = | 
holds in good approximation and thus (^) =1^ = 343. 
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Fig. 6. Resolution study of the maxima in the central shock col- 
lision of model z4.24A^d. 
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Fig. 7. One-dimensional resolution study of the maxima in the 
central shock collision of model s4. lOA. 



by another increase by about a factor of 145. Thus a maximum 
compression of about 2,900 can be achieved and in the case of 
■y = |, relevant in our case, those numbers are even bigger 

In our simulations the shock that propagates inwards into the 
WD due to the He detonation is followed by an extended wave of 
matter that is also flowing inwards. Therefore the case might lie 
somewhere in between these two idealized situations: The shock 
is amplified by the geometrical focusing eff'ect and thus com- 
presses a significant part of the core that is then burned at higher 
densities. On the other hand the matter succeeding the in-flowing 
shock leads to even higher densities after the shock is reflected. 
Taking all those properties into account, it becomes clear why 
the ignition conditions for a core detonation can be fulfilled so 
well in our simulations. Provided that the shock is symmetric 
enough also an off'-center ignition at a point with much smaller 
initial density like in model z4.24A_2d^56 is possible in this 
scenario. 

As the minimum shock surface that we can reach depends 
on the grid resolution, the maximum temperature and density 
values that are achieved through the geometric amplification in 
the shock implosion naturally also do so. In our simulations a 
successful second detonation ignition thus depends on the reso- 
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Table 4. Results of three complete double detonation simula- 
tions. All nucleosynthetic abundances are given in solar masses. 
Mioss is the mass loss over the grid boundaries. 
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Fig. 8. Time evolution of the mass fractions of model 
z4.24A_2d_128. 



lution and it seems likely that all the models in Table [3] would 
explode provided that they were well enough resolved. 

Therefore, the results given in Table |3]have to be interpreted 
as a conservative estimate of the possibility of a detonation ig- 
nition. If at a certain grid resolution the ignition conditions are 
fulfilled, a second core detonation must happen. However, if the 
ignition conditions are not fulfilled a core detonation can still not 
be excluded and might be achieved with increasing resolutiorj^ 

3.3. Three double-detonation supernova simulations 

As a first step to compare the simulations of sub-Chandrasekhar 
mass models with actual SN la observations, measurable quanti- 
ties like the nucleosynthetic abundances of an explosion have to 
be determined. To this end, two complete double detonations of 
model z4.24A until f = 2.0 s have been simulated. One of them 
was set up in 3D with 256^ cells for the full WD. For compari- 
son with previous results, the second, 2D rotationally symmetric 
simulation was also done with this low resolution (128 x 256 
cells). 

In order to test the influence of symmetry in the initial flame 
configuration in the He shell on the explosion abundances, be- 
sides these maximally asymmetric cases a 2D double detona- 
tion of the spherically symmetric model s4.10A was also simu- 
lated. The results of all three simulations are given in Table |4] 
For model z4.24A_2d_128 also the time evolution of the mass 
fractions of the species is plotted in Fig.[8]and will be discussed 
in the following. 

At first, the "^He in the shell is burnt completely to ""^Ni. But 
the mass of the latter is somewhat reduced due to some mass 
loss through the boundaries of the computational domain. This is 
clearly visible also in the decrease of the total mass. The reason 
for this is that the co-expanding grid was set up to track a mass 
shell at 0.85 Mq. Since the matter above this shell expands with a 



' Note that in this case the minimum detonator volume may be the 
limiting constraint. 



larger velocity, some of it leaves the numerical grid. The tracking 
of a bigger mass shell would reduce this mass loss, but it would 
also reduce the resolution of the colliding shocks in the core. 

After about 1 .2 s the collision of the shock waves resulting 
from the He detonation gives rise to conditions for a detonation 
in the C + O core. In the corresponding grid zones, this detona- 
tion was artificially triggered by setting up a second level set. At 
the high densities near the center the C/O mixture is burnt com- 
pletely to ^^Ni in NSE with a small amount of "^He. But the "^He 
recombines to ^^Ni after a short time due to the density decrease 
caused by expansion. Farther away from the center at smaller 
densities the burning product is ^**Si and some low-density mat- 
ter at the border of the core remains unburnt. In this way 0.40 Mq 
of ^''Ni and 0.51 Mq of^^S i are produced after about 1 .5 s and 
0.01 Mq of C/O remain unburnt. 

The results of the 3D case z4.24A_3d_128 are very similar to 
the 2D one. However, 0.05 Mq more ^^Ni were produced. The 
main reason for this difference is the loss of different amounts of 
mass over the boundaries of the computational domain, which 
amounted to 0.08 Mq in 2D and only 0.05 Mq in 3D. The expla- 
nation for the smaller mass loss is a simple geometric effect: In 
3D cartesian coordinates there is more volume outside the star 
that lies on the numerical grid than in 2D cylinder coordinates. 
For example, the largest distance from the center is the length of 
the space diagonal ( V3 a, when a is the side length), whereas in 
2D it is the length of the diagonal in the rz-plane ( 's/la). 

Due to the spherical symmetry in the model s4.10A_2dq_128 
He burning and shock convergence are faster. Therefore all the 
reactions finish already after about one second, but according to 
Table |4] almost the same nucleosynthetic abundances are pro- 
duced. This indicates that the density distribution in the core in 
case of shocks converging asymmetric to the density distribution 
differs only slightly from the perfect spherical symmetric case. 
As the now following core detonation propagates supersonic, the 
final abundances are independent of the position of the ignition 
point. 

In all three cases considered the ^*Ni-mass and with it the 
explosion energy are rather large. In fact, they reach the range 
of (weaker) normal SNe la. However, these numbers have to be 
taken with care since the density limit of 3 ■ 10^ gcm"^ adopted 
here for separating burning of C/O to ^^Ni and ^**Si, respectively, 
is quite uncertain, but will have an impact on the final abun- 
dances. 

We could not calculate synthetic spectra or light curves of 
our simulations yet. First order spectra and light curves are 
given, e.g., in Livne & Arnett ( 1995| l for models very sim- 
ilar to our z4.24A case. For comparison with their results 
the distributions of the species of models z4.24A_2d_128 and 
s4. 10A_2dq_128 in radial velocity space (averaged over all an- 
gles) are presented in Fig. [9] Both are very similar except for a 
minor difference in the unburnt C/O, which is more abundant in 
the spherically symmetric case. The velocity ranges are compa- 
rable to model 6 of |Livne & Arnett] ([T995), which had a similar 
setup (Mc/o = 0.8 Mq, Mhe = 0.17 Mq), but resulted in a larger 
''^Ni mass (0.648 Mq). Also, the maximum velocity of ^^Ni and 
the minimum velocity of ^**Si from core burning are shifted to 
significantly larger values. This could be caused by a difference 
in the transition density to NSE that might also be responsible 
for the difference in ''^Ni mass. 

In order to test possible asphericities of the ejecta in model 
z4.24A_2d_128 the abundances are plotted separately for the z-, 
r-, and -z-directions in Fig. 
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Despite the off-center explosion, 
the differences in the velocity distribution are minor Only in -z- 
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Fig. 9. Mass fractions of the species (averaged over all angles) of 
models z4.24A_2d_128 (solid line) and s4.10A_2dq_128 (dashed 
line) in radial velocity space at f = 10 s. 
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direction a shift in the maximum velocities of the ''^Ni from the 
core detonation and of the ^'^Si is visible. 



4. Summary and conclusions 

We have studied the sub-Chandrasekhar model of SNe la 
by means of a series of two-dimensional and a few three- 
dimensional simulations with different initial conditions. The 
numerical scheme used is based on the PPM algorithm, and 
the propagation of the detonation front is modeled applying 
the level set technique. This novel implementation allowed for 
an accurate treatment of the hydrodynamic features of the sub- 
Chandrasekhar model (such as shock waves and thin detonation 
fronts) in multiple dimensions. 

With our generic l-M© model we performed a parameter 
study involving different He masses and ignition geometries of 
the initial He detonation. We find that a detonation in the He 
shell can clearly trigger a second detonation in the core and at 
least for the mass configurations studied, the double detonation 
seems to be a very robust process, which works without any 
"fine-tuning" of our model parameters. In almost all of the sim- 
ulations performed, the ignition conditions for a core detonation 
were reached at or near the center of the WD as a result of the 
convergence of the shock from the He shell detonation and in 
the few other cases the failure is likely a result of the finite nu- 



merical resolution. The high maximum densities and tempera- 
tures that are needed for detonation ignition and the fact that the 
shock is not diluted too much when propagating inwards against 
the pressure gradient are made possible through a geometrical 
shock amplification effect that appears in spherical symmetry: 
As the surface of the shock decreases its strength has to increase. 
The high compression that can be achieved theoretically in this 
way, especially if full convergence is reached and the shock is 
reflected (cf. Sect. |3.2[ ), could even allow successful double det- 
onations with considerably smaller He masses. The shock am- 
plification that is reached, however, turns out to be resolution 
dependent on the numerical grid as it is coupled to the small- 
est resolvable shock surface. Taking this into account it is very 
likely that all our models would explode, if they were simulated 
with sufficiently high resolution. 

The question of whether an incineration could also happen 
directly at the core-shell interface, was not addressed here. This 
could of course prevent the spherical shock convergence mecha- 
nism from playing a role, at least in parts of the parameter space. 
Therefore, edge-lit detonations are postponed to a separate study. 

The complete double detonation simulations that were per- 
formed for the models z4.24A and s4.10A resulted in ^^Ni- 
masses of about 0.40 to 0.45 M^. Are the studied events thus 
good candidates for normal SNe la? Most likely not. Starting 
with a He-shell detonation, all of our models show a signifi- 
cant amount of rapidly expanding ^''Ni in the outer layers. In the 
observed spectra of normal and sub-luminous events, however, 
this has never been observed (cf. Branch et al.| ( |1982] l; |Branch| 
(| 1984| i ;|Woosley et al.|fi986|l;|Amett| ( |1997[ ) and also flie discus- 
sion in L ivne & Arnett( ( 1995| l). The only exception is the super- 



luminous SN 1991T. Thus our solar-mass models most likely 
will not be able to explain normal or sub-luminous SNe la and, 
given the robustness of the explosion mechanism of the model, 
the only conceivable explanation for the lack of observations of 
corresponding SN la events is that the progenitors are not real- 
ized in nature - or are very rare. 

For a better agreement with observations a reduction of the 
He-shell mass would be the most obvious choice. However, then 
a significantly larger core mass would most likely be required, 
because otherwise the He would not detonate. In this case the 
core density would be much higher resulting in a very large Ni 
mass. This would again not be a candidate for normal SNe la, 
but it might be a promising candidate for objects like SN 1991T. 

Further work could cover model improvements like a more 
realistic treatment of nuclear reactions including the calculation 
of real reaction rates depending on the actual thermodynamic 
state of the burnt matter This would make an investigation of the 
onset and the explosion dynamics of the core detonation possi- 
ble. Also, more extended parameter studies, especially towards 
the lower end of possible He masses, would be interesting. For a 
comparison with observations the calculation of synthetic light 
curves and spectra would also be desirable. 
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